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The mode dependence of current-induced magnetic excitations in spin valves is studied theoret- 
ically. The torque exerted on the magnetization by transverse spin currents as well as the Gilbert 
damping constant are found to depend strongly on the wave length of the excitation (spin wave) . An- 
alytic expressions are presented for the critical currents that excite a selected spin wave. The onset 
of macrospin (zero wavelength) vs. finite wavelength instabilities depends on the device parameters 
and the current direction, in agreement with recent experimental findings. 



PACS numbers: 75.70.Cn,72.25Mk 
I. INTRODUCTION 

Less than a decade ago Bergeni and Slonczewski^ ar- 
gued that an electric current sent through multilayers of 
normal metals (N) and ferromagnets (F) can excite the 
ferromagnetic order parameter and even reverse the mag- 
netization. The theoretical predictions have been con- 
firmed by many experiments on F|N|F nanostructured 
spin valves. '^'^ The physics of collective ferromagnetic ex- 
citations driven by out-of-equilibrium quasi-particles is 
complex and fascinating. In magnetic memories current- 
induced magnetization switching might turn out to be 
superior to its magnetic field driven counterpart. 

Current-induced magnetic excitations are driven by 
the spin transfer torque acting on the magnetic order 
parameter when a spin current polarized normal to the 
magnetization is absorbed by the ferromagnetii*SiSi& The 
transverse spin current extinction is a quantum mechan- 
ical dephasing effect between electrons at the Fermi en- 
ergy with different precession lengths. In Co, Ni and Fe 
this happens on an atomistic length scale. By conser- 
vation of angular momentum the absorbed spin current 
acts as a torque on the ferromagnetic condensate. At a 
critical spin current, this torque becomes strong enough 
to set the magnetization into motion, possibly leading to 
a complete reversal of the magnetization direction. 

One discussion that started with the prediction of the 
current-induced magnetization dynamics remains to be 
settled. Berger defines the critical current at the onset 
of spin wave excitations, whereas Slonczewski considers a 
rigid coherent rotation of the whole magnet ( "macrospin 
model). In the latter scenario, the critical current corre- 
sponds to a torque that just overcomes the Gilbert damp- 
ing. In our view, the differences in these pictures are to 
some extent semantic, since the macrospin model is iden- 
tical to the lowest energy spin wave. The physical ques- 
tion addressed here is the wavelength of the spin wave 
that is most easily excited. We find that there is no uni- 
versal answer and that the preferential excitation mode 
is a complicated function of device parameters and cur- 
rent direction. Nevertheless, our theory agrees well with 



experiments that observe both type of excitations >^ 

Most theories of spin transfer torques and critical cur- 
rents are based on macrospin precessions in spin valves^*^ 
Indeed, sufficiently small magnetic clusters support a sin- 
gle domain magnetization and the magnetic field induced 
magnetization reversal is well described by a coherent ro- 
tation according to the Stoner-Wohlfarth modelii From 
a theoretical point of view, the macrospin models valid- 
ity is essential for understanding the non-linear physics 
underlying the entire magnetization dynamics by dynam- 
ical systems theory and the probabilistic treatment in the 
"presence" of noise. ^ In a single ferromagnetic film sand- 
wiched by normal layers, the torque on the macrospin 
domain necessarily vanishes. However, biased N|F|N 
structures are unstable with respect to spin waves with 
shorter wavelengths In spin valves, we may there- 
fore expect a competition between macrospin and shorter 
wave length spin waves excitations. Experiments on spin 
valvesi^ have indeed been interpreted in terms of both 
types of excitations, depending on the current direction 
(and the spin accumulation pattern). It is our purpose 
to understand and model these data in order to assess 
the dependence of the excitation modes on device pa- 
rameters. Another motivation to study the competition 
between different excitation modes is the need to find 
criteria for the breakdown of simple models for the mag- 
netization that provide a guide for the necessity of full- 
fledged micromagnetic calcul ations !^ 

Our work extends Refs. Ildlll on single layers to 
spin valves. We derive analytical expressions for bias- 
driven spin transfer torques, enhanced Gilbert damping 
constants^ and critical currents for magnetic excitations 
as a function of wave vector. We predict a rich "phase 
diagram" in the current vs. magnetic field plane. 

Our paper is organized as follows. Section ^ intro- 
duces our model for the spin valve pillar. We also present 
the assumptions used in computing the charge current, 
the spin current, the spin transfer torques and the en- 
hanced Gilbert damping governed by spin pumping. Fur- 
thermore, a schematic description of the necessary in- 
gredients for the calculations is outlined. Section IIIII 
presents the computational details for readers interested 
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FIG. 1; Spin-valve a)with circuit resistance elements sketched 
in b). 



in a deeper knowledge on how our results have been ob- 
tained. Our results are presented and discussed in sec- 
tion IIVI In that section, we also compare our results 
with the measurements by the NYU/IBM collaboration 
and find a semi-quantitative agreement. We conclude our 
paper in section 



II. MODEL AND ASSUMPTIONS 

We consider N|F|N|F|N spin-valve pillars in Fig. ^ 
in the semiclassical transport regime where magneto- 
electronic circuit theory applies.''^ Disregarding spin flips 
in F is allowed for sufliciently thin F films and consider- 
ably simplifies the calculations. We also consider normal 
metals thinner than the spin-diffusion length. A possi- 
bly spatially dependent pillar cross section, not discussed 
below any further, can be treated by simply scaling the 
resistance parameters. 

We compute the spin transfer torques, the enhanced 
Gilbert damping and the critical current for magnetic 
excitations in the following way. The charge and spin 
accumulation pattern and the charge and spin cur- 
rent through the system is computed by treating the 
N|F|N|F|N spin valve layer by layer. Our calculations 
includes scattering due to bulk impurities within each 
layer and scattering due to the band-structure mismatch 
between adjacent metals, e.g. interface resistances. Both 
bulk impurity scattering and interface mismatch are im- 
portant in e.g. understanding the giant magnetoresis- 
tance effect in magnetic multilayers. Bulk impurity scat- 
tering is taken into account by solving the diffusion equa- 
tions in the bulk of each material. Spin and charge cur- 
rent through four interfaces between the normal metals 
and the ferromagnets are determined by the interface 
scattering matrix between the adjacent materials which 
gives rise to the interface resistances that we take into 
account. 

We will show below that our results for the spin trans- 
fer torques, the enhanced Gilbert damping and the crit- 
ical current for magnetic excitations can be expressed in 



terms of the resistance parameters of the electronic cir- 
cuit in Fig. Let us first introduce the resistances 
that arc important for the calculations and the results: 
Starting from the left in Fig.^ we introduce a bulk resis- 
tance per unit area for the bulk of the left normal metal, 
~ PNitNn where p^i is the resistivity of the normal 
metal and tjq^ is the thickness of the normal metal layer 
in the transport direction. Similarly, starting from the 
right, we introduce the bulk resistance per unit area for 
the bulk of the right normal metal, Rr = PN^tN^ a-nd the 
bulk resistance per unit area for the bulk of the middle 
normal metal, Rm = PN^tN^. The interesting physics 
arises due to the spin-dependent scattering within the 
bulk of the ferromagnetic layers and at the interfaces be- 
tween the normal metals and the ferromagnetic layers. 
The spin-dependent resistances Rf {i = 1,2, s =T, J.) in 
the ferromagnetic elements consist of interface and bulk 
contributions: i?| = 2Rp^j^ + Qptpi, where Rp.^j^ is the 
spin-dependent resistance per unit area of a single F|N 
interface assumed identical on both sides of film i. gp. 
and tp. are its spin-dependent resistivity and thickness. 
When the magnetizations are coUinear the spin depen- 
dent resistances in the ferromagnets with the resistances 
in the normal metals are sufficient to describe the trans- 
port of charge and spin in the systems. When the mag- 
netizations are non-coUinear, the flow of spins directed 
perpendicular to the magnetization in the ferromagnets 
must also be quantifled. The mixing conductance G]^ de- 
termines the absorption of spins transverse to the mag- 
netization direction and consequently the spin transfer 
torque^ We find it convenient to introduce a mixing re- 
sistance R]^ = l/(2Gj^) so that all our results can be 
presented in terms of resistances. 

The charge and spin currents are not only induced 
by the applied bias, but also by spin pumping due to 
a moving magnetization that opposes the magnetization 
dynamics^ Effectively, when a layer magnetization pre- 
cesses, it acts as a source of a spin current transverse 
to the magnetization direction and proportional to the 
precessional frequency, which should be included in the 
self-consistent computation for the charge and spin cur- 
rents in magnetoelectronic circuits. When spin- waves are 
excited in one of the layers, the sources of spins are non- 
uniform. This leads to a wave-vector dependence of the 
enhanced Gilbert damping. These contributions can be 
computed by using the same equivalent circuit. Fig. 
as for the bias voltage induced charge and spin current 
and also depend on the same resistance parameters of 
the circuit. The wave- vector dependent enhanced Gilbert 
damping turns out to be governed by the magnitude of 
the spin-pumping current out of the ferromagnet and the 
resistance between the ferromagnet and either the volt- 
age reservoirs or other spin sinks for the emitted spin 
current. 

We wish to compute the critical currents for the mag- 
netic instabilities in driven spin-valve pillars by lineariz- 
ing the equations that govern and the magnetization dy- 
namics close to the equilibrium configurations. In this 
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way, we can find the critical currents for magnetic exci- 
tations, but we cannot unambiguously predict what hap- 
pens when the current is increased beyond the critical 
current. To the latter end, one needs to go beyond the 
linear instability regime and study the full micromag- 
netic behavior of each individual ferromagnet as well as 
the transport of spin and charge self-consistently. This 
is beyond the scope of the present work. We would 
like to emphasize that even in the linear regime the dy- 
namics of the ferromagnets and the transport of spin 
and charge have to be treated self-consistently. It is 
the self-consistent feedback of the spin-flow from a non- 
uniform domain back into the ferromagnet that destabi- 
lizes the monodomain configuration in a the single layer 
ferromagnet 

We are mostly interested in spin- valves that have a rel- 
atively large in-plane magnetic field along the magneti- 
zation direction as e.g. studied by the NYU/IBM collab- 
oration. Therefore, we proceed from equilibrium configu- 
rations of monodomain magnetizations in parallel or anti- 
parallel to each other. Out of equilibrium, spin waves 
can be excited transverse or parallel to the transport di- 
rection. For sufficiently thin ferromagnets, the critical 
current for the transverse excitations are lower than the 
critical current for longitudinal excitationsiiSiii Conse- 
quently, we consider small transverse instabilities of the 
magnetization direction in Fourier space, Sm.j^{q) and 
i5m^(q), where q is a two dimensional vector in Fourier 
space transverse to the transport direction and compute 
the perturbed charge and spin current through the sys- 
tem. The perturbed spin currents induce spin transfer 
torques on the ferromagnetic layers. Also, a precessing 
ferromagnet emits spin current into the adjacent materi- 
als than can enhance the Gilbert damping. We compute 
whether the magnetic instabilities are stable or unsta- 
ble and find the critical current for magnetic excitations. 
Longitudinal magnetic excitations along the transport di- 
rection (as studied e.g. in Ref. 11) arc thus disregarded. 
We will comment on this thin-layer assumption when ap- 
plying our theory to the NYU/IBM experiments in Sec. 

El 



transverse direction, mi = mi(p) and m2 = m2(p). 
The charge current along the N|F interface thus depends 
on the transverse coordinate p. The charge currents on 
the two sides of the first and second ferromagnet-normal 
metal interface are 
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are the charge voltage drop over the ferromagnets 1 and 



2. The spin voltage drop A\f' = vJ'J^ 



'^2m defined analogously. Charge 
current is conserved on traversing the ferromagnet, which 



Vjf and 



is reflected in the conditions = and jj! 

jjji''. The spin-dependent conductances Gf (i = 1,2 
and s =t, i) consist of interface and bulk contributions: 
1/G| = Rf = 2Rpjj^ + gp.tpi, where Rp./^ is the spin- 
dependent resistance per unit area of a single F|N inter- 
face assumed identical on both sides of film i. g^, and tp. 
are its spin dependent resistivity and thickness. Similarly 
to the charge current, we express the spin current density 
incident from the left (middle) normal metal and going 
into ferromagnet 1 as ji; (jim) and the spin current den- 
sity incident from the middle (right) normal metal and 
going into ferromagnet 2 as j2m (j2r)- The spin-current 
density is expressed in terms of parallel (||) and trans- 
verse components (±), 
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III. CALCULATION OF CURRENT 



The parallel components are 



We introduce charge V^'^'' and spin V'^''^ potentials in 
the normal metals close to the N[F interfaces. These po- 
tentials depend on the two-dimensional coordinate p = 
(x, y) along the interface transverse to the transport di- 
rection z. In our notation, ju (jim) denotes the charge 
current density incident from the left (middle) normal 
metal and going through ferromagnet 1 with magnetiza- 
tion along unit vector mi, and j2m (j2r) denotes simi- 
larly the charge current density incident from the middle 
(right) normal metal and going through ferromagnet 2 
with magnetization along unit vector m2. The direc- 
tions of the magnetizations also depend on the trans- 
verse coordinate p allowing spin-wave excitations in the 
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The component of the spin-current parallel to the mag- 
netization direction is conserved on traversing the ferro- 
magnet which is reflected in the conditions jjj'"'' — Jim'' 

and j2m^ = i2r"'' ■ The transverse components of the spin- 
current densities are not conserved on traversing the in- 
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where it is assumed that the mixing conductance is the 
same for the interface between the left (middle) normal 
metal and ferromagnet 1 (2) as for the interface between 
the middle (right) normal metal and ferromagnet 1 (2), 
and the small imaginary part of the mixing conductance 
G]^ has been disregardedi^ The spin-torque densities on 
ferromagnets 1 and 2 are determined by the absorption 
of the transverse spin currents: 
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^From the expressions for the spin current densities JSJ, 
©, O and (jHJ and the flow of spins transverse to the 
magnetizations (|13|) . H14II . (|15|l and Hlt)|) . we can express 
the spin transfer torque densities in terms of the spin 
accumulations in the normal metals: 
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We consider normal metals thinner than the spin- 
diffusion length, but spin-flip may not be neglected. Eqs. 
O, ©, ®, and gjl and ©, ©, O and © should 
be matched to solutions of the diffusion equations for 
charges and spins 



(21) 



where qnt] is the resistivity and Irj the spin-diffusion 
length in each normal metal layer with index rj — l,m,r. 
We denote the direction along the spin valve pillar, e.g. 
the transport direction, z. The charge and spin currents 
in the transport direction are 
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The thicknesses of the left, middle and right normal met- 
als are L/, Lm, and Lr, respectively. We assume that the 
left and right normal metals are attached to reservoirs 
where charges and spins are in local equilibrium. We 
choose local coordinates, zi, Zm, and z^, for each normal 
metal so that 1) in the left normal metal = corre- 
sponds to the contact between the left side normal metal 
and the left reservoir and zi = Li corresponds to the 
Ni|Fi interface, 2) in the middle normal metal z„i ~ 
corresponds to the right side of the Fi|Nm interface and 
Zm = corresponds to the left of the Nm|F2 interface. 



and 3) in the right normal metal Zr = corresponds to 
the right side of the F2|Nr interface and Zr = Lr corre- 
sponds to the boundary between the right normal metal 
and the reservoir. 

Conservation of spin and charge requires continuity of 
spin and charge flow in the bulk of the normal metals 
close to the left and to the right of the N|F interfaces 
and through the interfaces: 
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On the outer, left and right, normal metal layer, the 
boundary conditions are fixed charge voltage and zero 
spin potentials defining the magnetically active device 
region. At the left reservoir, we have 



Vi^'\zi ^ 0) ^ V/2 , V\'\zi =0) = 
and at the right reservoir, we have 

Vy\zr - Lr) - -V/2 , Vi'\zr = L,) = . 



(27) 



(28) 



The charge and spin accumulations are uniquely deter- 
mined by Eqs. C3l-®i l|21|) . H22|l . conservation of charges 
and spins close to the N|F interfaces in the normal metal, 
and the boundary conditions on the outer normal metal 
layers. 

We start from a parallel or antiparallel spin valve con- 
figuration in which magnetizations and spin accumula- 
tions are coUinear to x. In this initial configuration, 
the transport equations are easy to solve since both 
the spin accumulations and magnetizations directions are 
collinear throughout the whole circuit, and e.g. the trans- 
verse components of the spin currents vanish. The cur- 
rent through the system then follows similarly to Ohm's 
law for resistances. 

We now consider small instabilities in the magneti- 
zation normal to the x axis in Fourier space: m^ — 
x-l-^m^(q). To linear order in the excitations, the charge 
current and voltage become j^'^^ — j^'^"-' -I- (Sj*-'^-'(q) and 
l/(c) = yico) _|_ sycf^f^y y^^Q decompose the linearized 
spin potentials into longitudinal and transverse parts, 

expressing the spin transfer 

torque density n = jjf ^ -i^^ (r, = - j^f ^) in 
terms of contributions from electrons hitting Fi from the 
left and right: ti = tu + Tim {t2 = T2m + T2r): 
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Solving the diffusion equations in each normal metal, we 
can also find the charge currents to linear order in the 
instability close to the N|F interfaces 



-^j(?)(0) = G, 



tanh qLi 
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cosh qL„ 
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(33) 
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(35) 
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where l/Gi = Ri ^ PnM, l/G^ = i?^ = p^^L^, 
and 1/Gr = Rr ~ PN^Lr are the resistance of the left, 
middle and right normal metal layers, respectively. The 
expressions for the spin-currents are similar by replac- 
ing j'"^) — > j'"^ and F^^' V^'*). The equations for the 
desired charge and spin potentials are can be obtained 
by noticing that expressions for 5j'^'^'^ and Jj^'*"-' form a 
closed set of equations independent of 5va.±_ with trivial 
solution: SV'^""' = JV^^'ID = and Sj^'"'^ = Jjt^H' = 0. 
Only transverse spin currents j^**"*"^ are therefore induced 
by small magnetization fluctuations. Finding the solu- 
tions of the linear equations for the unknown 5'V\i , 



(si) 



(5V, 



2m 



and ^V2^ is straightforward, but te- 



dious. We present the results in the next section. 



IV. RESULTS AND DISCUSSION 

In order to manage the complex analytical spin transfer 
torque expressions, we introduce the following standard 
notation for the ith fcrromagnet: average resistance 



R* = iR]+Ri)/4 



and polarization 



P, - 



Rj + Rl 



(37) 



(38) 



The resistance contrast between parallel (P1P2 > 0) and 
anti-parallel (P1P2 < 0) configurations is 



AR = i?^P -RP = 



miR*2 IP1P2I 



{Ri + RI + R,n + R*2 + Rr 
The polarization of the current is 

p ^ RlPi + R*^P2 

Rl + R'^ + Rjn + R2 + Rr 

The charge current is denoted 



(39) 



(40) 



A. Spin transfer torques 

The spin transfer torque on ferromagnet 1 (2) has con- 
tributions from electrons hitting it from the left (middle) 
normal metal and from the middle (right) normal metal, 

Tl = Til + Tim {t2 = T2m + T2r)- Tli = CuSmj^j^''^ On 

the first ferromagnet exerted by electrons coming from 
the left and T2r = >C2r^ni^j^'^) on the second ferromag- 
net due to electrons from the right is given by: 



Cii 



Rl im-Mq )] 

RiMq)+Rl 



n 



Pt. 



(41) 



Expression for C2r are obtained by changing sign and 
substituting / by r and 1 by 2. Here f{x) = tanh(x)/a;, 
h{q) = f[xi{q)\, and Xi{q) = [q^ + I^^Y/Hni, where 
t^i is the left normal metal thickness and its spin- 
diffusion length. In particular /;(0) « 1 (and simi- 
larly for fm, fr, and Xm, Xr)- These torques acting 
on a single ferromagnetic layer vanish in the long wave- 
length limit as in Ref. ITol The new physics is con- 
tained in the wave vector-dependent torques Tim and 
T2m- In the expressions Tim = (£im(5m5'--|-/Cim(5m^)j(^' 
and T2m = ('C2mi5m^ -I- lC2m5T^i)j^'^\ thc spin transfer 
torques due to spin currents between the two ferromag- 
nets are found as: 



t^lr 
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Rl + R^RVxI + Rm {rV + RV) /fn 
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where, 

711 - -RUPiRI + P2R*2) , (43) 

n2 = RV[PiR*i{R*2 + Rm + Rr) - P2R*2{Rl + R*l)] , 

(44) 
(45) 



713 = niRl^ /Rm + n^Rm/R, 

introducing R* — Ri + RI + Rm + R2 + Rr- The expres- 
sions for £2™ is similar to that for Cim, with an overall 
sign change and the substitution I r and 1 2. K,im 
and /C2m govern the dynamic coupling between the ferro- 
magnets, but do not affect the instabilities of individual 
ferromagnets and are therefore not discussed here. In the 
limit of long-wavelength excitations, g — s- 0, our results 
reduce to previous ones.^'-'^'* 

Let us discuss the spin transfer torques H41() and H42() 
in simple limits. 

(i) When the ferromagnet 2 and the right normal metal 

are absent in the circuit, e.g. when i?2 = Oi -^2 = 
and Rr = 0, we naturally obtain a similar result as the 
spin-torque on the first ferromagnet exerted by electrons 
coming from thc left (|41|> ^ 
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RrnifmiO) - .fmj q)] 
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(46) 



When, additionally, the system is symmetric around fer- 
romagnet 1, P; = Rm, the total spin transfer torque on 
ferromagnet 1 vanishes, Cu + Cim = OiiS 
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(ii) For symmetric junctions in a parallel magnetic con 

figuration Ri = R^, R\^ = RI^ , Rl = R*2, and Pi ^2- ^ enhances the Gilbert damping ai = a\ 
For macrospin {q — 0) excitations 



ferromagnet i into neighboring normal metals cx Cj^mj x 
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and a2 = a 
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(47) 
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illustrating that the spin transfer torque is exerted by 
electrons coming from ferromagnet 2 that hit ferromag- 
net 1 from the right and is governed by the polarization 
of the entire spin-coherent region Pt (Ref. M). At short 
wavelengths, 



^11 — Pt — tT J 
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so that 



Cl — Cil + Cir 
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(49) 
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For the dynamics of a single ferromagnet, it is the asym- 
metry in the diffusion process to the left and to the right 
of the ferromagnet that determines the sign of the short 
wavelength spin transfer torqueiiS We see that for spin- 
valves, it is the asymmetry of the complete device that 
determines when short wavelength modes cannot be ex- 
cited. When Rm = 2Ri, small wavelength modes cannot 
be excited. This condition {Rm = 2Ri) corresponds to 
a completely symmetric spin valve, where two identical 
effective "ferromagnets" consists of e.g. i?/, ferromagnet 
1 and half of ferromagnet Rm and half of ferromagnet 
Rm, ferromagnet 2 and Rr- In the limit Ri = there is 
no short wavelength spin transfer torque from electrons 
that hit ferromagnet 1 from the left and vice versa when 
Rm = we find Tun{q 0) ^ 0. On the other hand, the 
macrospin, torque in symmetric junctions only depends 
on the global polarization for symmetric junctions^. We 
thus provided proof of the conjecture in Ref. [T^ that the 
spin transfer torque Cim can change its sign as a function 
of q when Rm > 2Ri for symmetric spin-valves. In that 
case, the magnetization moves as a macrospin for one 
current direction but short wavelength spin waves are 
excited when the current is reversed. In general, we find 
that when the normal-metal bulk resistance asymmetry 
is smaller than the interface spin-mixing resistance the 
macrospin spin transfer torque is larger than the short 
wavelength spin transfer torque. For strongly asymmet- 
ric structures, Cim is negative also for long- wavelength 
excitations 



B. Enhanced Gilbert damping 

The critical current for the onset of magnetic insta- 
bilities depends also on spin-pumping by a moving mag- 
netization that opposes the dynamics.^ Spin pumping by 



damping parameters in the isolated ferromagnets. Spin 
pumping via the left normal metal to the left reservoir 
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where 7* is the gyromagnetic ratio. Mi and Vi are the 
magnetization and volume of ferromagnet i, and Rk = 
h/e^ is a conductance quantum. The enhancement due 
to spin pumping via the middle normal metal to ferro- 
magnet 2 comprises a new result: 
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and a'lm can be understood in terms of an ef- 
fective resistance against pumping spins out of the fer- 
romagnet. Spin pumping to the left is limited by the 
wavelength-dependent effective total resistance R\^ + 
Rifi- In the long- wavelength limit the resistors are in 
series, rI^ + Ri. In the short-wavclcngth limit, the effec- 
tive resistance is reduced due to the inhomogeneous spin 
distribution. Spin pumping to the right, in the long- 



wavelength limit governed by R[^ + Rm + RI^ is also 
reduced in the short- wavelength limit. For the second 



ferromagnet, we compute a2 = a^-^"^ + a'^i + a'2 
identical expression as above with 1 <-> 2 and I <-> 



with 



C. Critical Currents for magnetic excitations 

Including spin transfer torques and enhanced damping, 
the magnetization dynamics obey a generalized Landau- 
Lifshitz- Gilbert equation: 
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= -jnii X Heff - 
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The effective magnetic field Hofr is the functional deriva- 
tive of the total magnetic energy U due to external 
magnetic field, magnetic anisotropy and spin-wave stiff- 
ness. To lowest order in the excitation amplitude U — 
-mf^ ■yiH + Ki{ui-5mlY + K2{n2-5ml f, where the 
applied magnetic field H is assumed to be along the ini- 
tial magnetization direction. Ki{q) and K2{q) are wave 
number dependent spin-wave stiffnesses that for g = 
reduce to the anisotropy constants along axes Ui and U2. 
A configuration in which ferromagnet 1 is aligned par- 
allel with the external magnetic field becomes unstable 
when j^'^' > , whereas an antiparallel magnetization 
becomes unstable when j^'^^ < . We compute the crit- 
ical currents at T = requiring that the damping torque 



is exactly canceled by the spin transfer torque: 



JO) 



^'ll + a'lr 



2e „ cj± 
Mitpi—, 

n 7* 



(55) 



introducing the ferromagnetic resonance frequencies 
Lj^ = H ± {Ki + K2)/2. and Cf^ denote spin trans- 
fer torques computed for positive/negative polarization, 
Pi = ± |Pi|. When the torques change sign as a function 
of wave vector the ferromagnet can be unstable against 
macrospin excitations when the current is flowing in one 
direction, but spin-waves are excited when the current is 
reversed, as pointed out above. Similar expressions hold 
for j2 the critical current for ferromagnet 2, j^. 
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FIG. 2; Phase diagram for ferromagnet 1 and 2, starting in a 
parallel configuration, as a function of current. 
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D. Comparison with NYU/IBM experiments 

We are now in position to discuss the experiments in 
Ref. ^3 A good measure of the efficiency of the spin 
transfer torques is the critical current dji/dH at high 
magnetic fields when anisotropics Ki and K2 do not 
play a role.^^ The magnetoresistance, enhanced Gilbert 
damping, spin transfer torques and switching current 
all depend strongly on the detailed device parameters. 
We adopt resistance parameters collected by the Michi- 
gan State University group for Co/Cu systems, i.e., 
Qq^ — 75 nilm, bulk polarization /3 = 0.46, ^^^u = 
6 nfim, Pco/Cu ~ ^-^^ fOm~^, and an interface po- 
larization 7 — 0.77. The mixing conductance that 
agree with first-principles band-structures calculations 
and ferromagnetic resonance experiments is R\^^/q^ = 
1.0 ffim"^ (Ref. The magnetically active region is 
10Cu|3Co|10Cu|12Co|35 Cu, where the numbers denote 
the thickness of the layers in nanometers. The right nor- 
mal metal thickness is chosen smaller than the geometri- 
cal since the nano-pillar widensii^ 

Before going into details, we summarize our results 
for the NYU/IBM samples in Fig.|21in semi-quantitative 
agreement with the experiments. Starting in a parallel 
configuration there is a macrospin instability for ferro- 
magnet 1 at positive currents. The macrospin excitation 
leads to an antiparallel configuration for larger currents,"* 
which should be confirmed by spatially dependent micro- 
magnetic calculations. A further increase in the current 
leads to a spin- wave instability in ferromagnet 2. Sim- 
ilarly, for negative current, we predict a spin- wave in- 
stability for ferromagnet 2 (Ref. IT^ . The theoretical 
Gilbert damping constant is plotted in Fig. O for a bulk 
Gilbert damping a[?^ = 0.003 for both ferromagnets. 
Spin-pumping is important, giving rise to a strongly en- 
hanced Gilbert damping constant of the thin ferromagnet 
that increases with wave vector. Both spin-pumping and 
spin transfer torque increase with the interface to volume 
ratio. Ferromagnet 1 is easier to excite but its enhanced 
damping partially compensates this effect, thus allowing 
excitations of both ferromagnets to compete. Numerical 
results for the critical current dji/dH at high magnetic 
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FIG. 3: Upper: Gilbert damping in thin ferromagnet 1 (up- 
per curve), thick ferromagnet 2, and in bulk (dashed line) 
vs. qtMm- Lower: Critical current dj/dH in units of 
10*'^Am~^T~^ of thin ferromagnet 1 (solid line) and thick 
ferromagnet 2 (dashed line) vs. qtMm- 



fields are displayed in Fig. Ofor the thin ferromagnet 1 
and the thick ferromagnet 2 for the parallel configuration. 
For positive currents, experiments find switching of the 
thin ferromagnet presumably to an antiparallel configu- 
ration when dj'-''^/dH = 0.6 x lO^^ Am^^T^i (Ref. 0). 
This value is extracted from the line in their contour plot 
of (PV/cPl showing the onset of magnetic excitations. 
Indeed, for positive currents macrospin excitations oc- 
cur first with lowest critical current that agrees well with 
experiments. Furthermore, we expect spin-wave excita- 
tions for the thick ferromagnet 2 at opposite currents, 
again in good agreement with experiments.-^ Addition- 
ally (but not shown here) we predict that after the thin 
ferromagnet 1 switches, ferromagnet 2 also becomes un- 
stable against spin-wave excitations for positive currents, 
supported by experiments as well. 

Residual quantitative discrepancies between experi- 
ments and theory are believed to be caused by uncer- 
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tainties in the material and devices parameter and details 
in the micromagnetic structure. We have used measure- 
ments for e.g. the bulk resistivities from the MSU group 
and it is likely that due to a different sample preparation 
technique, the NYU/IBM samples have different resis- 
tivities. Furthermore, the thickest ferromagnet is f2nm 
in the NY/IBM measurements, which is at the boundary 
where spin- wave excitations along the transport direction 
become of importance in this layer, which reduces the ac- 
curacy of the predicted value of the critical current in the 
thick layer»ii Given these uncertainties, we consider the 
semi-quantitative agreement encouraging. 

V. CONCLUSIONS 

In conclusion, we report a theory of macrospin vs. 
spin-wave excitations in spin valves that explains recent 
observations using only independently determined ma- 
terial and device parameters. The rich phase space of 
magnetic excitations is classified in terms of macrospin 



and finite wavelength spin-wave excitations that depend 
on the resistance distribution in the magnetically active 
region. For symmetric junctions, macrospin instabilities 
are strongly favored. Finite wavelength spin-wave exci- 
tations are pronounced in asymmetric spin-valves with 
relatively high normal metal resistances comparable to 
that of the ferromagnets. Since the results are in agree- 
ment with recent experiments, we are confident that our 
insights should be helpful to explore the magnetization 
dynamics in the full parameter space spanned by cur- 
rents, external magnetic fields, and device design. 
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